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We propose a simple discrete model to study the nonequilibrium fluctuations of 
two locally coupled 1+1 dimensional systems (interfaces). Measuring numerically 
the tilt-dependent velocity we construct a set of stochastic continuum equations 
describing the fluctuations in the model. The scaling predicted by the equations 
are studied analytically using dynamic renormalization group and compared with 
simulation results. When the coupling is symmetric, the well known KPZ exponents 
are recovered. If one of the systems is fluctuating independently, an increase in the 
roughness exponent is observed for the other one. 
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The dynamics of equilibrium and nonequilibrium interfaces has attracted much atten- 
tion recently The continuum theory proposed by Kardar, Parisi and Zhang (KPZ) |2| 
provided a very successful analytic approach for different problems, ranging from spin sys- 
tems to growth models. However usually the dynamics of interest is dominated by a strong 
coupling fixed point, making the perturbative approach inconclusive. In such a situation 
simple discrete models allow us to obtain detailed information on the dynamic scaling. 

In this paper a generalized single step model is introduced to study the nonequilibrium 
dynamics of two coupled 1 + 1 dimensional systems. Investigating the tilt dependence of 
the growth velocity we find that the set of nonlinear equations describing the fluctuations 
of the two surfaces ho(x,t) and hi(x,t) is 

d t h = u dlh + \ (d x h ) 2 + ^d^d^x + ^(d^) 2 + r] (la) 

d t hx = Pid 2 x hi + Ait^/ii) 2 + jAhod^ + <pi(d x h ) 2 + ^ (lb) 

where d t and d x denote the partial derivatives with respect to t and x. The noise is assumed 
to have zero mean and uncorrelated {rji(x,t)rji(x',t')) = Di5(t — t')5(x — x') with % — 0, 1. 
The scaling of the fluctuations are determined numerically from the discrete model and the 
results are compared with the predictions of the dynamic renormalization group (DRG) 
calculations. A slightly modified version of the model allows us to study the evolution of 
a growing interface perturbed by a nonequilibrium field. The increase of the roughness 
exponent of the surface observed numerically is supported by the analytic predictions. 

Equation (1) represents the most general form describing the fluctuations of two coupled 
nonequilibrium systems with local dynamics. It describes the evolution of two interfaces 
which are moving with the same mean velocity in an inhomogeneous media and are coupled 
by some local interaction. More generally it describes the growth of an interface ho(x,t) in 
the presence of a nonequilibrium field h\(x,t), assuming that the interface and the field are 
strongly affecting each others behavior. In this respect it might provide a suitable analytic 
framework with which to describe the recent surface growth experiments in which the control 
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of other physical factors was not possible. Such a perturbing field might be the density or 
pressure of the fluid during fluid displacement in porous media || , the concentration of the 
nutrient in bacterial growth [[| , or the density of the mesoscopic particles and impurities in 
the imbibition experiments ||. 

On the other hand the study of equation (1) is of general interest in the context of 
the recent efforts to understand the general properties of the nonequilibrium stochastic 
systems. Recently Ertas and Kardar (EK), studying the motion of a single flux line in 
random environment, have presented detailed investigation of (1) for 70 = Ai = if\ = 
0. Imposing an additional condition, (po = 0, the problem reduces to the convection of a 
passive scalar (T = hi(x,t)) in a random velocity field (v = Vh Q ) 0. For 7! = ipi = 
equation (1) describes the fluctuations of an interface h (x,t) evolving under the influence 
of an independent nonequilibrium field hi(x,t). || 

The Model - Our goal is to introduce a simple model to study the scaling of the two 
coupled interfaces. For this we use the mapping between the single-step model and the 
driven hard-core lattice-gas on a one- dimensional chain |J. In the mapping each step along 
the surface ho(x, 1) is associated with a site on the chain H Q . A site j is said to be occupied 
(i.e. H (j) = 1) if we have an upward step on the interface h (x = j,t), and is empty 
(H (j) = 0) otherwise. Two one-dimensional lattices, H (j) and Hi(j) with (j = l,L), 
are filled with probability |. First choose randomly a site j on one of the lattices (for 
example H ). If it is filled (i.e. H (j) = 1), the quantity $ = [H (j — 1) + H (j) + H (j + 
1) — 1.5] [Hi(j — 1) + Hi(j) + Hi(j + 1) — 1.5] determines the direction in which the chosen 
particle wants to move. If $ > the Ho(j) particle will try to move to the left, otherwise to 
the right. The step can be completed only if the neighboring site in the chosen direction is 
empty. Then choose randomly an another particle on the other lattice and repeat the same 
procedure. The unit of time is defined as L trials per lattice, where L is the system size. 
Periodic boundary conditions are imposed on both of the lattices. 

The quantity [H (j — 1) + H (j) + H (j + 1)] is related to the local slope of the interface 
ho(x = j, t), thus with a positive $ the two interfaces have the same slope and is generating 
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a growth on the site j, while a negative $ corresponds to different slopes and results in a 
decrease of the height. 

We are interested in the scaling of the fluctuations characterized by the dynamic scaling 
of the width M 



w 



t 



(L, t) =< (h t (x, t) - h t (t)) 2 >= L 2 ^f( — ) (2) 



L z 



where x% is the roughness exponent for the interface hi(x,t), and the dynamic exponent 
z% describes the scaling of the relaxation times with length; hit) is the mean height of the 
interface at the moment t and the < ... > symbols denote ensemble and space average. The 
scaling function fi has the properties fi(u — > 0) ~ u 2l3i and fiiu — > oo) ~ consti, where 

Tilt- dependence - In order to identify the relevant nonlinear terms determining the scaling 
behavior in the model we have measured the tilt dependent velocity ||10|| , imposing a global 



slope on the interfaces h (x, t) and hi(x, t). A positive slope is induced on h (x, t) by increas- 
ing the concentration of the H particles, a negative slope is a result of a decrease in concen- 
tration. Fig. 1 shows the tilt-dependent velocity v and v i, where Uj = J2j=i,L(dthi(x = j, t)). 
The velocities which may explain the observed behavior and are compatible with the 
x — > —x symmetry of the model have the form: v o = k — ao{d x hod x hi) — bo(d x ho) 2 and 
Vi — k + ai(d x liQd x hi) — bi(d x h ) 2 , where k, a iy bi are positive constants. The tilt in h (x, t) 
induces a negative mixing term (d x hod x hi) in t> and a positive one in v\. But since in the 
model the influence of the Hq{x, t) to hx(x, t) is identical to the influence of hi(x, t) to ho(x, t) 
(symmetric coupling), without tilt one expects the same sign for both mixing terms. The 
symmetric coupling is responsible for the presence of a <^o(<9x^i) 2 term in vq and a \i(d x hi) 2 
term in i>i, which can be observed if the tilt is in hx(x,t) instead of h (x,t). Based on 
these measurements we conclude that equation (1) contains all the relevant nonlinear terms 
determining the scaling of the two interfaces, and it describes the proposed model if Aj < 
and tpi < 0. Higher order terms might be responsible for the unusual behavior observed for 
large tilts, but they are in fact irrelevant concerning the scaling. 



Dynamic Renormalization Group Analysis - The scaling behavior of the equation (1) in 
general can be investigated using one- loop DRG |f7,|H|. Rescaling the parameters x — > e e x, 



t — > e £z t, and hi — ► e iXi hi, we obtain the following flow equations for the coefficients: 

dvg _ r .K"i Ag-Dp Ti^ogi , i^i7o , ^PiDq JqDi , 

d£~ U0[Z Z+ u [ vl + 2*? ) + 2v (v + Vi) [ u + v x ] 

(^o + v x ) 2 vq v x 



^^_ 2 , _ 1 + W + ^£A ] + m (3) 

d£ z^o z/ z/i(z/ + z^i) v{ 



dX . r 2<pi-D 7oP , ix,7oA - ^i , 2<p <Pi 7o7i 
Ao[^ + Xo - 2] + ■ ( )(- 



dt Vq + Vi Vq V\ Vq Vq + V X V X V + ^1 

^7o r . m , -^l w 2 <Pi D o ToAw (7o-2Ai) A y? . 7i<^0x . 
-77 = 7o * + Xi - 2 + ( (7o + -= — ) + 

dl Vq + Vx V Vx Vq + Vx Vq ZVx 

2(p Dx 7iA)w 7i-2A A 7o </Wi 
7o : + ~z 

Vx Vq V + Vx 2V Vx 

d^Po ro , ,71 A) 2^ £>iw 7i-2A 70 - 2Ai 
— = 2xi - Xo + ^ - 2 + : Wo + 7o : 

dl V + Vx V Vx Vq Vq + V\ 

The flow equations for the other five missing coefficients can be obtained from (|3]) replacing 
vq <-> vx, Dq <-> Di, X <r-> Ai, 7 <-> 71, y2 <-> and Xo ^ Xi- 

Obtaining the exponents from (^) is not straigforward because of the large number of 
parameters involved, but important results can be obtained by making use of the nonper- 
turbative properties of (1) combined with the direct integration of (|]). 

If 2lpqVqDx = JiViDq and 2ipxVxD = ^qVqDx (fluctuation-dissipation (FD) subspace) the 
joint probability f| 



p[h (x,t),hx(x,t)\ = exp(- J dx[^-(d x h ) + 7p^( d M D 

is a solution of the Fokker-Planck equation following from (1). This provides us the exact 
exponents Xi — 1/2- The nonlinear terms do not renormalize, resulting in the scaling relation 
z + Xi = 2. Direct integration of (§) shows that for positive (pi and 7^ the flow converges 



to the FD subspace (see Fig 2), thus resulting in the super diffusive exponents Xi — 1/2 and 
2 = 3/2. A change of variables ho — > —ho and hi — > —hi indicate the same exponents if all 
the coefficients of the nonlinear terms are negative. 

As the velocity measurements indicated, Aj < and <pi < in the model. If 7$ is 
also negative, the FD subspace dominates the behavior. The DRG does not provide exact 
results in the 7« > case; direct integration of (H) indicates a strong coupling fixed point, 
with diverging A,, ji and ipi. 

Fig. 3 shows the scaling of the width w 2 (L, t) ~ t 2 ^. The exponents were determined 
using saturated systems of size L = 50, 100,200,400 and unsaturated systems of size L = 
1000,2000,4000. The roughness exponent Xi was obtained from the best collapse for the 
time dependent width (0). All the simulations indicated super diffusive exponents, giving 
as result x% = 0-52 ± 0.03 and $ = 0.32 ± 0.02. Concluding this section we note that the 
simulations are in perfect agreement with the predictions of the DRG for 73 < 0, so it is very 
likely that this is the sign of the mixing term in the model. But since no analytical results 
are available for 7, > 0, we can not rule out the possibility that for this sign the fluctuations 
of the two systems are also characterized by the superdiffusive exponents. 

In the following two sections we study a variant of the proposed model, in which one of 
the interfaces fluctuates independently. This corresponds to the above mentioned situation, 
when an interface ho(x,t) grows in the presence of an equilibrium or nonequilibrium field 
hi(x, t). The fluctuations of the field is assumed to be independent of the growing interface. 
The presence of the field results in an enhancement of the roughness of the growing surface. 

The A = <Po = <fi = 71 = case - The hi interface fluctuates independently of h 
in the following version of the model: The randomly chosen particle on the Hq chain will 
try to move in the direction determined by the sign of but the one chosen on H 1 will 
try to move left with a probability p and right with probability 1 — p, independently of 

On Hi we have the single step model with evaporation. If p 7^ 1/2, the fluctuations of 
the hi interface are described by the KPZ equation, thus 71 = tpi = in (1). For p = 1 
measurements on the tilt dependence of vq indicate the absence of the Ao and ipo terms: the 
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velocity does not change if we impose a tilt either on ho or on hi, the only influence is the 
decrease of V\ for a tilt on h 2 , as a result of the nonvanishing negative Ai term. But these 
results do not indicate the absence of the 70 term: the fact that it has no influence on the 
velocity may come from the vanishing contribution of < d x hod x h\ > , which is acceptable 
considering that hi fluctuates independently of ho. The exponents determined from the 
model are xo — 0.64 ± 0.03 and (3q = 0.33 ± 0.02, together with the known exponents of 
the KPZ equation: xi = 1/2 and (3\ = 1/3 (see Fig 4). These exponents are in good 
agreement with those obtained by EK from the direct integration of the equations (1) for 
Ai < and 70 > [ 13| . The flow equations (|3|) show that Ai scales to zero, thus the DRG 



is not conclusive for these coefficients. It indicates only that z > Zi, which agrees with 
the numerical findings. A support for the positivity of 70 comes from the DRG result for 
7o < 0, the one-loop exponents being xo = 3/4 and j3 = 1/2, considerably larger than those 



observed numerically in the model [14 



The Ao = Ai = y?o = V^i = 7i = case ~ If i n the second variant of the model we 
choose p = 1/2, the hi interface fluctuates in equilibrium with Ai = tpi — 71, leading to 
the Edwards- Wilkinson (EW) exponents Xi — 1/2 and Pi = 1/4 The situation in the 
model is the following: the hi interface fluctuating in equilibrium influence the fluctuation 
of the ho interface through the mixing term (d x hod x hi). The velocity measurements again 
confirm the absence of the (p^ Aj terms. The simulations indicate the nontrivial exponents 
Xo = 0.68 ±0.02 and (3 = 0.34 ±0.02, together with the EW exponents for h x . The DRG in 
the present form can not be applied because of the different scaling of the time in the two 
interfaces (i.e. z 7^ Zi). 

Discussion - In principle the existence of a term (d^hi) in (la) (and (c^o) i n (lb)) 
can not be excluded neither from symmetry considerations nor from velocity measurements. 
These terms can be transformed away at the expense of a cross-correlated noise in (1) (i.e. 
the new noise has the form r] = Co + and ?7i = Co _ a Ci) where a is a constant and Co 
and Ci are uncorrelated noises). Further study is necessary to determine the influence of a 
such a noise on the scaling and on the DRG flow. 



The results presented below are just a small subset of the possible behavior following from 
(1) and (^). Further investigations are necessary to understand the scaling for other choices of 
the coefficients. In this paper we have focused mostly on those parameter values which were 
accessible through the studied discrete models. Determining the relevant nonlinear terms 
using tilt-dependent velocity measurements, we have shown that the model is described by 
two coupled nonlinear Langevin equations of form (1). For symmetric coupling the model 
allowed us to determine the scaling exponents which were consistent with the exact results 
obtained from the DRG. The study of two particular cases were also possible using a modified 
version of the model. In the first one a nonequilibrium field ho was influenced by hi, whose 
oscillations are described by the KPZ equations, while in the second case hi was fluctuating 
in equilibrium. Both cases leaded to considerable enhancement of the roughness exponent 
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FIGURES 

FIG. 1. The tilt dependent velocity vq (empty symbols) and v\ (filled symbols) for sys- 
tem sizes L = 125 (circles), L = 250 (triangles) and L = 500 (squares). Each set of data 
points was obtained by averaging over 250,100 and 25 independent runs. The inset shows the 
simplest polynomial in the tilt (d x ho) of the form vq = 0.031 — 0.34(d x ho) — 0.12(d x ho) 2 and 
v\ = 0.031 + 0.02(921/10) — 0.02(0 x /io) 2 , which may account for the observed behavior in the vicinity 
of zero. Higher order terms may be responsible for the complicated behavior observed for large 
tilt. 

FIG. 2. The convergence of the flow (3) to the fluctuation-dissipation subspace for positive 
coefficients, starting from different initial conditions. For simplicity the symmetric coupling is 
presented, i.e. uq = v\, Dq = Di, Ao = Ai, 70 = 71, <^o = fi and Xo = Xi- This choice is motivated 
by symmetry between ho(x,t) and h±(x,t) in the model. Direct integration of (3) shows that the 
symmetric subspace is also stable for small perturbations. 

FIG. 3. Symmetric coupling: The scaling of the width with time for systems of size 
L = 50, 100, 200 and 400. An average over 50000, 10000, 5000, and 1500 independent runs were 
taken. The main figure shows the scaling after an intrinsic width [12] of magnitude 0.225 was 
extracted from the data. The slope of the straight part gives Po = Pi = 0.32 ± 0.01. The inset (a) 
shows the same data without extracting the intrinsic width. The roughness exponent xo = 0.52 
was determined from the best collapse of the data, shown in inset (b). 
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FIG. 4. The Ao = (fi = 71 = case: The scaling of the width with time for systems of size 
L = 200, 400 and 800. An average over 1000 independent runs were taken. The main figure shows 
the scaling after an intrinsic width [12] of magnitude 0.225 was extracted from the data. The slope 
of the straight part gives Po = Pi = 0.33 ± 0.02. The roughness exponent x was determined from 
the best collapse of the data, shown in inset (a). Inset (b) shows the scaling of the height-height 
correlation function < \hi(x,t) — hi(x + l,t)\ 2 >~ l 2x ' for the interfaces ho(x,t) (empty symbols) 
and h\(x,t) (filled symbols), for the same system sizes as in the main picture. The slope of the 
straight line is 2xo = 1-28. 
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